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ABSTRACT 

We study the way Lyot coronagraphs with unapodized entrance pupils respond to small, low order 
phase aberrations. This study is applicable to ground-based adaptive optics coronagraphs operating 
at 90% and higher Strehl ratios, as well as to some space-based coronagraphs with intrinsically higher 
Strehl ratio imaging. We utilize a second order expansion of the monochromatic point-spread function 
(written as a power spectrum of a power series in the phase aberration over clear aperture) to derive 
analytical expressions for the response of a 'band-limited' Lyot coronagraph (BLC) to small, low 
order, phase aberrations. The BLC possesses a focal plane mask with an occulting spot whose opacity 
profile is a spatially band-limited function rather than a hard-edged, opaque disk. The BLC is, to 
first order, insensitive to tilt and astigmatism. Undersizing the stop in the re-imaged pupil plane 
(the Lyot plane) following the focal plane mask can alleviate second order effects of astigmatism, at 
the expense of system throughput and angular resolution. The optimal degree of such undersizing 
depends on individual instrument designs and goals. Our analytical work engenders physical insight, 
and complements existing numerical work on this subject. Our methods can be extended to treat the 
passage of higher order aberrations through band-limited Lyot coronagraphs, by using our polynomial 
decomposition or an analogous Fourier approach. 

Subject headings: instrumentation: adaptive optics — space vehicles: instruments — techniques: 
image processing — astrobiology — circumstellar matter — planetary systems 



1. INTRODUCTION 

Lyot coronagraphy ijLvotJ Il930l Il939j) has enjoyed a 
resurgence because interest in discovering and char- 
acterizing extrasolar planets has been stimulated by 
advances in ground-based adaptive optics (AO), as well 
as support from space agencies for extrasolar plane- 
tary science missions. There is currently one extreme 
adaptive optics (ExAO) coronagraph being used to 
conduct a complete complete surve y for faint compan- 
ions a nd disks around nearby st a rs IIOpd enhcimer^t^dj 
120031 12001 IDigbv et all l200l IMakidon et alJ 120051) . 
New coronagraph designs to enable the detection 
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and characterization o f extra sol ar planets a bound 
(INisenson fc Papaliolio'sl 120011: lAime et alJ I200S : 
Kuchner fe Traubl 12002] iSoummer et alJ l2003albl: 
Aime fc Soummerl 12003 iKasdin et all 120031 iSoummeJ 
2005). Many of these designs would suppress light 
from the central star around which a planet orbits, 
as long as the telescope and instrument optics and 
stops were perfect, and simple Fourier optics theory 
was accurate enough to predict instrument behavior 
in these high contrast regimes. Less-than-ideal image 
quality is a stumbling block which impedes attainment 
of this science goal. Kasdin (private communication) 
has demonstrated in the laboratory that Fourier optics 
appears to be valid down to contrast ratios of 10~ 7 . 
We study the way Lyot coronagraphs dedicated to 
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imaging extrasolar Jovian and terrestrial pla nets respond 
to sm all, low order phase aberrations. iPerrin et alJ 
( 2003) demonstrated that for ground-based adaptive op- 
tics systems delivering Strehl ratios of at least 90%, a 
second-order expansion of the imaging system's point- 
spread function (PS F) response to monochromatic light 
from a point source ijSivaramakrishnan et al.ll2002[) was 
sufficient to model the effects of phase aberrations on 
the PSF. We analyze the case of unapodized telescopes 
in high dynamic range regimes, as is relevant to the prob- 
lem of tolerancing optics on space-based coronagraphic 
telescopes or calibrating wavefront sensing non-common 
path errors on ground-based ExAO instruments dedi- 
c ated to finding and characte ri zing e xtrasolar planets. 

iLlovd fc Sivaramakrishnanl l|2005f) derived the re- 
sponse of a simple band-limited coronagraph to small 
tilt errors, using an analytical method combined with 
an expansion of the pupil plane field strength (de- 
fined later) in terms of the phase aberration over the 
pupil. That work produced insight into the way tilt 
error causes light to leak through a coronagraph that 
is designed to suppress all light from a perfectly flat, 
on- axis wavefront. Such an analytical treatment de- 
velops a qualitative understanding of coronagraphic re- 
sponse to phase errors. Insights developed there apply 
to classical Lyot coronagraphy. We extend the methods 
of ILlovd fc Sivaramakrishnanl ((2005) to understand how 
low order aberrations propagate through a band-l imited 
Lyot coronagraph. While earl ier stu dies (e.g.. | Malbetl 
Il99fit ISivaramakrishnan et alJ l200ll iGreen et alJ 12003ft 
explored some aspects of coronagraphy on imperfect 
wavefronts, they did not provide analytical tools required 
to understand some of the optical problems of designing 
a coronagraph to achieve the 10 7 photometric dynamic 
range needed by a space mission dedicated to discover- 
ing and characterizing extrasolar Jovian planets, or the 
10 10 contrast ratio that e xtrasolar terrestrial planet dis- 
cover y may well demand l)Breckinridge fc Qppenheimerl 
l2004j) . Our analytical approach is a start along the route 
to understanding coronag r aphic data describe d in e.g., 
lOrroenheimer et alJ l|2004l) : iDigbv et all l|2004fl . 

A band-limited coronagraph is a perfect coronagraph 
design in that simple Fourier optics modelling suggests 
that it will prevent all incoming, on-axis light from 
reaching the final coronagraphic focal plane. We define 
light leak as the fraction of incident on-axis energy that 
reaches the final coronagraphic focal plane. Our concern 
is restricted to diffractive light leak rather than light re- 
distributed within the instrument by scattering off im- 
perfect optics. We also concentrate on small aberrations 
of the wavefront that forms the first image. For instance, 
we treat tilts that are small compared to the angular res- 
olution of the system. Such tilts decenter the star behind 
the focal plane mask's occulting spot without moving it 
out from behind the spot — which is itself at least a 
few resolution elements in diameter. Such aberrations 
in the wavefront prior to the focal plane mask result in 
light leaking through a supposedly perfect coronagraph, 
reaching the final coronagraphic image plane, thereby re- 
ducing the dynamic range of the instrument. We find 
that small amounts of tilt and astigmatism do not in- 
crease the light leak through the coronagraph, whereas 
the same amounts of coma and spherical aberration do. 
However, detrimental second order effects of astigmatism 



can be reduced by undersizing the Lyot stop. Our anal- 
ysis shows that defocus and tilt produce characteristic 
intensity distributions in the Lyot plane, which suggests 
ways this signal can be used to align a high dynamic 
range coronagraph, and calibrate the offset between the 
focal plane mask location in the final coronagraphic im- 
age plane, thereby improving the coronagraph's astro- 
metric accuracy. Similar analyses on the various corona- 
graph designs being developed for the purpose of detect- 
ing Jovian and Earth-like planets around nearby stars 
will contribute to an informed selection of a design, in 
addition to providing insight into how to improve these 
designs. 

2. BASIC CORONAGRAPHIC THEORY 

A simple mathematical de scription of the fundamen- 
tals of_JjVot_corona^ra£hy llLvotl ll939T) can be found 
in ISivaramakrishnan et all (120011) lAime fc Soummerl 
( 2002)7 lLlovd fc Sivaramakrishnanl 1)20051) . and refer- 
ences therein. 

2.1. Monochromatic coronagraphic theory 

Here we briefly recapitulate our basic monochro- 
matic Fourier optics formalism, wh ich fol - 
lows that of |LiQ^xLfc_SiYar a,makrishnanl |2005); 
Sivaramakrish nan fc Llovdl (l2005fl . A more d etailed 
treatment can be found in iBorn fc WoM ©93). We 
recollect that a plane monochromatic wave travelling in 
the z direction in a homogenous medium without loss of 
energy can be characterized by a complex amplitude E 
representing the transverse {e.g., electric) field strength 
of the wave. The full spatio-temporal expression for the 
field strength is £: e ( lKZ -"*) ; where lu/k — c, the speed of 
the wave. We do not use the term field to denote image 
planes — the traditional optics usage — we always 
use the term to denote electromagnetic fields or scalar 
simplifications of them. The wavelength of the wave 
is A = 2tt/k. The time-averaged intensity of a wave 
at a point is proportional to EE*, where the average 
is taken over one period, T = 2ir/u;, of the harmonic 
wave. The phase of the complex number E represents 
a phase difference from the reference phase associated 
with the wave. A real, positive E corresponds to an 
electric field oscillating in phase with our reference 
wave. A purely imaginary positive value of E indicates 
that the electric field lags by a quarter cycle from the 
reference travelling wave. Transmission through passive, 
linear filters such as apertures, apodizers, and so forth, 
is represented by multiplying the field strength by the 
transmission of these objects which modify the wave. 
Again, such multiplicative modification is accomplished 
using complex numbers to represent phase changes 
forced on the wave incident on such objects. 

We assume that scalar Fourier optics describes ou r 
imaging system (cf. iBreckinridge fc QpDenheimerl l2004'l : 
image field strengths are given by the Fourier transform 
of aperture (or pupil — we use the two terms interchange- 
ably) illumination functions, and vice versa. 

A telescope aperture is described by a transmission 
function pattern A(x), where x = (x,y) is the location 
in the aperture, in units of the wavelength of the light 
(see Fig. The corresponding aperture illumination 
describing the electric field strength in the pupil (in re- 
sponse to an unaberrated, unit field strength, monochro- 
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Fig. 1. — The essential planes and stops in a coronagraph. The 
entrance aperture is A, the direct image at B falls on a focal plane 
mask whose transmission function is m(k). The re-imaged pupil 
plane D, after being modified by passage through a Lyot stop with 
a transmission function N(x), is sent to the coronagraphic image 
at F. A, D, and E are pupil planes, and B, C, and F are image 
planes. 
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Fig. 2. — A sketch of the classical Lyot coronagraph (left), and 
the 1 — sine "top-hat" band-limited coronagraph operating on a 
perfectly flat incoming wavefront. The top row shows the same 
aperture function used in both examples, viz., a clear unobscurcd 
entrance aperture. The second row shows the field strength at the 
image plane. The third row shows the field strength immediately 
after passage through the focal plane mask. The classical Lyot, 
on the left, has a hard-edged, opaque mask. This band-limited 
coronagraph has a mask that is opaque only at its center. The 
fourth row shows the field strength in the Lyot pupil plane before 
a Lyot stop is applied, and the final row shows the field after the 
Lyot stop in the Lyot pupil plane. The band-limited coronagraph 
blocks all incoming on-axis light if the wavefront is unaberrated. 



matic incident wave) is Ea — A(pc). From this point 
onwards we drop the common factor Ee^ lfiZ ~ ult ^ when 
describing fields. The aperture intensities (EaE a ) for 
two coronagraphic designs are shown in Fig. 2 (top row). 
The field strength in the image plane, Eb — a(k), is the 
Fourier transform of Ea, where k = (k x , k y ) is the image 
plane coordinate in radians. Because of the Fourier rela- 
tionship between pupil and image fields, k is also a spatial 
frequency vector for a given wavelength of light. We refer 
to this complex-valued field a as the 'amplitude-spread 
function' (ASF), by analogy with the PSF of an optical 



system. The PSF is aa* . Our convention is to change the 
case of a function to indicate its Fourier transform. We 
multiply the image field Eb by a mask function m(k) to 
model the focal plane mask of a coronagraph. The image 
field immediately after this mask is Ec — m(k) Eb- The 
electric field in the re-imaged pupil after the focal plane 
mask is En, which is the Fourier transform of Ec- We 
use the fact that the transform of the image plane field 
Eb is just the aperture illumination function Ea itself, 
so the Lyot pupil field is Ed — M(x)*Ea, where * is the 
convolution operator. 

If the Lyot pupil stop transmission is AT(x), the electric 
field after the Lyot stop is Ee — N(x)Ed- The transform 
of this last expression is the final coronagraphic image 
field strength: Ep = n(k)*[m(k) Eb]- 

The aperture illumination function with phase aberra- 
tions is 



A(x)e' 



(x) 



(1) 



In this paper we look at the way small phase aberrations 
0(x) pass through a band-limited coronagraph. Our ap- 
proach is to expand the exponential in equation Q in 
powers of </>(x) , so the quantity that determines the rate 
of convergence of the expansion is the largest excursion 
of 4> from its mean value over the clear aperture. 

2.2. The Lyot plane field 

A coronagraphic stop is modelled by multiplying the 
image field by a focal plane mask transmission function 
m(k), so the electric field strength just after the focal 
plane mask is a a f, er m. When describing a Lyot corona- 
graph it is sometimes helpful to introduce a 'mask shape 
function' w by the definition 

m(k) = 1 - w(k). (2) 

When w is unity at the origin the mask is opaque at its 
center. We can write the electric field at the Lyot plane 
as 

E Ly ot(x) = A aber *(S(x.) - W(x)) 



(3) 



Understanding the subtleties of this equation in the case 
of high Strchl ratio imaging is essential to understanding 
how small phase aberrations cause light to leak through 
a coronagraph designed to produce perfect on-axis image 
suppression with completely unaberrated on-axis light. 

2.3. Band-limited coronagraphs 

IKuchner fc Traubl 1I2OOI use a focal plane mask shape 
function w which is band-limited. This means that there 
is a minimum positive value of b such that the mask 
function's FT, W, satisfies the property 



W{x) = if |x| > b. 



(4) 



The bandwidth or bandpass of w is b (note that b is 
actually a physical distance in pupil space). We se- 
lect the mask function in its transform (pupil) space, 
(x,y), rather than in physical (image) space, (/cj;,/^^), 
even though it is applied in the image plane in any real 
coronagraph. If the telescope diameter is D, then the 
characteristic scale of the mask function is D/b resolu- 
tion elements (a resolution element is X/D radians, A 
being the wavelength of the monochromatic light). This 
results in a focal plane mask about D/b Airy rings in 
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2.4. The top-hat band-limited coronagraph 

The simplest choice for W(x) is the 'top-hat' function 
U(x/d,y/d)/d 2 , where 

IL(x,y) = 1 for jorj < 1/2, }y| < 1/2, f . 

II(x, y) = elsewhere ^ ' 

(see Fig- EJ - The bandwidth of the corresponding mask 
function is d/2, A normalizing factor of 1/d 2 is applied to 
ensure that the area under W(x) is unity, thus ensuring 
that w(0,0) = 1; the focal plane mask is opaque at its 
center. A top-hat W produces a multiplicative image 
stop function 

m(k) = 1 — smc(dk x ) smc(dk y ). (6) 

Our interest such a choice is primarily didactic: this 
coronagraph leaks aberrated light in a manner similar to 
that of the more popular 'sawtooth' band-limited coron- 
agraph, viz., one with a mask function FT of 

W(x, y) = A(x/d, y/d)/d 2 = A(x/d)A(y/d)/d 2 , (7) 

where 

A(a?) = for \x\ < 1, , , 

A(a;) = elsewhere. ^ ' 

The sawtooth coronagraph has a focal plane mask de- 
scribed by m(k) = 1 — sine 2 (dk x ) sine 2 (dk y ). 

3. EXPANSION USING THE PHASE FUNCTION 

The phasor e 1 ^ can be expanded as 1 + i<\> — <f> 2 /2\ + ... 
for any finite value of the aberration <f>. Depending on 
the size of the aberration it is u seful to truncate this 
expansion at various orders {e . g.. iBloem hof et al 1l2001i 
iSivaramakrishnan et alJ l2002t iPerrin et al] 120031) . The 
first order perturbation of the field strength in the aper- 
ture is imaginary: i.e., the field due to phase aberrations 
is in phase-quadrature to the perfect aperture field. The 
second order aberration of the aperture field strength is 
in (anti-)phase with the perfect wave's field strength. 

3.1. Zernike polynomials and Cartesian expansions 

Zernike polynomials, which are a set of orthogonal 
polynomials on the unit radius 2-dimensional disk, are 
frequently used t o describe phase aberrations in opti- 
cal systems (e.g., iBorn fc Wollll^ lNolll97l . These 
polynomials are expressible in either circular coordinates 
(r, 9) or Cartesian coordinates (x, y). Expressing Zernike 
polynomials in Cartesian coordinates leads to a mathe- 
matical simplicity that promotes a better understanding 
of the physics of an imperfect band-limited Lyot coron- 
agraph. Because of this, we investigate coronagraphy on 
a square aperture in order to make integrals separable in 
x and y. 

We have the choice of of writing the phase as either a 
simple polynomial or an expansion in terms of Zernikes. 
In the latter case, we must choose the Zernike functions' 
normalized coordinates: we choose the disk on which we 
define our Zernike functions to possess an area equal to 
that of our square aperture. This choice results in extrap- 
olating the Zernike polynomials to a radial coordinate 
of 2/y/n ~ 1.128, which is beyond their usual range of 
validity. Doing so is a fruitful exercise if we restrict our- 
selves to low order aberrations. We emphasize that we 
do not invoke any of the Zernike polynomials' orthonor- 
mal properties to prove any result. Thus our apparently 



cavalier extrapolation of the polynomials beyond the do- 
main over which they are usually defined does not lead 
to grief. We only use them here because they are widely 
used in the optical literature. A more natural Fourier de- 
composition of the phase aberration over the aperture is 
of course possible, and may even preferrable to a Zernike 
decomposition for higher order aberrations. 

Our choice of the equal-area circle for our Zernike poly- 
nomials is motivated by the following considerations: if 
we use the circumscribing circle around our square aper- 
ture, we miss the steep increases of Zernike functions 
describing e.g., spherical aberration at the edge of the 
disk, which for the most part will lie outside our square 
aperture. If we use an inscribed circle, we include large 
areas in the corners of the apertures where the magni- 
tude of the third, fourth, and fifth order polynomials 
increase sharply, contributing features that we do not 
commonly associate with particular Zernike aberrations. 
Instead of either extreme, we stick to the middle ground, 
at the expense of mathematical purity: we miss some of 
the higher values of these Zernike polynomials where our 
circle extends outside the square aperture, but collect 
larger phase errors at the corners of our aperture, where 
the normalized Zernike arguments exceed unity. The 
net effect of this is that we can still talk about Zernike 
functions on our square aperture, in much the same way 
Zernike function aberrations on obscured apertures such 
as HST or Palomar are discussed, even t hough the annu- 
lar Zernike functions of iMahaianl l|1981l) should be used 
on annular apertures. 

4. PROPAGATION OF SIMPLE POLYNOMIAL PHASE 
ABERRATIONS 

At any location in the pupil plane, A a b er can be ex- 
panded in an absolutely convergent series in <f> for any 
finite value of the phase function: 

A aber = AA^ = A(l +i4>- 4> 2 /2 + ...). (9) 

If the phase over the aperture is expressed by 

oo 

^{x,y) = ^2,aiZi(x,y) 

l=1 oo (10) 

^ y CV-mn% y ; 
n— 0,m— 

where the Zt's are Zernike polynomials, and a%, a mn are 
constant coefficients, then equation @ produces first, 
second, and higher order terms in x and y. A measure of 
the magnitude of the phase aberration — e, the largest 
absolute value of the deviation of the phase <j) from its 
aperture-weighted mean — can be used to estimate the 
size of a particular order term in equation 10. Given any 
particular value of e one can estimate whe re the terms 
in th e above expansion become negligible l)Perrin et alJ 
12001 . For example, in the Advanced Camera for Sur- 
veys (ACS) coronagraph on the Hubble Space Telescope 
(HST), which sits in the aberrated beam, with about half 
a wave of spherical aberration at a wavelength of 800nm, 
one must use about seven terms of the above expansion 
to model the gross features of the aberrated field strength 
over the aperture. 
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4.1. Power leak through a top-hat band-limited 
coronagraph 

We can calculate how an arbitrary polynomial phase 
aberration of the form a nm x n y m propagates through the 
coronagraph to the Lyot plane, using this expansion. The 
calculation is straightforward for a square aperture be- 
cause the limits of integration and the basis functions are 
separable. From equation J2J we see that the Lyot field 
for such an aberration is described by 

ELyot(x,y) = A(x,y)a nm x n y m 

-TL{x/d, y/d)/d 2 *[A(x, y)a nm x n y m ] 

for the 'top- hat' band- limited coronagraph design. This 
calculation becomes difficult to perform analytically 
when the mask transform function has a more extended 
support, as in a classical Lyot design. 
We define the function 

L n , n (x,d) = Il{x/d)*x n /d, (12) 

which enables us to write equation 1)11(1 as 

EL V ot{x,y) = a nm x n y m 

L u ,n(x,d)Lu, m (y,d) K ' 

in the interior of the Lyot plane in the case of un- 
apodized apertures. We define the interior of the aper- 
ture as that part of the aperture (or image of the aper- 
ture, such as the Lyot plane) which is further than b 
from any point on the aperture boundary (where b is 
the bandwidth of the mask shape function). This is the 
area of the Lyot pupil which would not be obstructed 
by a Lyot stop in a perfect band-limited coronagraph 
(see iKu chner fc Traubll2002HLlo yd fc Sivaramakrishnanl 
120051 ISiva ramakrishnan & Llovd 120051 for further de- 
tail). The aperture function A(x) is unity everywhere in 
the interior of the aperture for an unapodized entrance 
pupil, and so does not appear explicitly in equation (|13|) . 

The 2.4 m HST ACS coronagraph has an image plane 
stop diameter of l'/8, or ~ 30 resolution elements at 
800 nm. Therefore the results we derive here for band- 
limited coronagraphs apply even though the ACS coro- 
nagraph is not band-limited, because its Lyot plane can 
be said to possess an 'interior' area which is further than 
a few -D/30 scale le ngths away from the annular aper- 
ture boundaries ( see ISivaramakrishnan et al.ll200l|, for a 
discussion of the natural scale length in the Lyot plane) . 

Explicit evaluation Lu in (x,d) is a matter of straight- 
forward integration: the first few values of n produce 

1 - L U fl(x,d) = 

x — Lu i(x, d) =0 

x 2 - Ln, 2 (x,d) = -d 2 /12 

x 3 -Ln3(x,d) = -d 2 x/A 

x 4 - L nA (x, d) = -d 2 x 2 /2 - d 4 /80 

x 5 -L n , 5 (x,d) = -5d 2 x 3 /6- d 4 x/16 

x 6 - L n6 (x, d) = -M 2 x 4 /A - M 4 x 2 /l& - d 6 /448. 

(14) 

The right hand sides of these equations are the first order 
residual Lyot field strength given a small phase aberra- 
tion of a n x n over the entrance pupil of a one-dimensional 
band-limited coronagraph. 

Low order phase aberrations, along with the resultant 
first and second order Lyot plane intensities are shown 



Tilt 


t GCU.3 




.□ 










1 o 




o 




! 


m 



Fig. 3. — Response of the 1 — sine 'top-hat' band-limited coron- 
agraph to tilt, astigmatism, focus, coma and spherical aberrations. 
Each aberration shown produces a 95% Strehl ratio direct image. 
Phase aberrations (top row) result in non-zero intensity in the Lyot 
pupil interior. Analytical first order (middle row) and second or- 
der (bottom row) approximations to these intensities are shown on 
a logarithmic grey scale, between 0.1 (white) and 10 -7 (black). 
First order power leak through the BLC due to tilt and astigma- 
tism are identically zero. Defocus, and second order tilt cause 
uniform illumination in the Lyot pupil interior, which produces a 
fainter version of the direct PSF (although with a larger resolution 
element due to Lyot stop undersizing) in the final coronagraphic 
image. 

in each column of Fig. [3J The band limit of this corona- 
graph is b = D/5. The fraction of incident power leaking 
through into the Lyot plane interior is presented in Ta- 
ble n Each power leak fraction is renormalized by divid- 
ing by the fractional coronagraphic throughput, which is 
the ratio of the areas of the optimal Lyot stop and the 
entrance ap erture. For a band limit of -D /5 this ratio is 
9/25 (e.g., Sivaramakrishnan et al. 2001). 

5. EXTENSION TO THE SAWTOOTH BAND-LIMITED 
CORONAGRAPH 

Designs similar to the sawtooth (sinc-squared focal 
plane mask function) coronagraph are being considered 
for the James Webb Space Telescope's Near Infrared 
Camera, as well as for a dedicated coronagraphic space 
telescope. In order to calculate how the sawtooth band- 
limited coronagraph leaks low order aberrations, we can 
define a function L\ q(x, d) using the analogy of equation 
03: 

L A<n (x,d) = A(x/d)*x n /d, (15) 

and carry out the same analysis as in section 14.11 We 
present the top-hat band-limited coronagraph's response 
to low order aberrations (Fig. |3J) ; the sawtooth corona- 
graph's response is qualitatively similar. The latter coro- 
nagraph's mask function's Fourier transform is a saw- 
tooth function, which also yields tractable algebraic ex- 
pressions for the Lyot pupil field strength for any order. 
For the sawtooth design the residual Lyot field in the in- 
terior of the Lyot plane has similar algebraic properties 
to that of the top-hat design. The powers of the relevant 
physical quantities for the 1 — sine 2 mask function of the 
sawtooth coronagraph are the same as those in equation 
1(14(1 . although the numerical coefficients differ slightly. 

6. SECONDARY MIRROR ALIGNMENT ERRORS 

When considering a simple on-axis two-mirror tele- 
scope, errors in positioning the secondary mirror result 
in aberrations described by a few low order Zernike poly- 
nomials. 
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Off-axis designs may have different sensitivities, but 
their phase aberrations due to secondary mirror tilt, 
despace, and decenter are still described by combinations 
of low order polynomials. 

6.1. Sensitivity to stop size 

The amount of energy leaking into the interior of the 
Lyot pupil is proportional to the fourth power of the 
bandwidth of the stop, or the inverse fourth power of 
the stop size expressed in resolution elements. This fol- 
lows from inspection of equation (|14|l . The leaked energy 
finds is way into the final image plane, unless the Lyot 
stop is adjusted to block some of this light. This strong 
dependence must play a significant role when choosing a 
coronagraph design given a desired contrast ratio derived 
from a scientific goal. 

6.2. Crosstalk between aberrations 

When two aberration, either simple monomials or 
Zernike polynomials, are present simultaneously, the 
Lyot pupil field contains energy due to a 'crosstalk' term 
in addition to effects from each individual term, even 
in the first order calculation. Cataloguing the various 
crosstalk energy leaks is beyond the scope of this work: 
we restrict ourselves to a few simple, traditional aberra- 
tions. 

6.3. The effects of tilt and defocus 

Equation (|14J) demonstrates that for a band-limited 
coronagraph there is complete cancellation of the central 
source even for small tilt errors (£n,i = 0). This particu- 
lar result has already been derived bvlKuchner fc Traubl 
i|2002ft : lLlovd & Sivaramakrishnan! (|2005^ : here it is de- 
rived as part of the band-limited coronagraph's response 
to a sequence of small polynomial phase errors. 

The way focus affects the field strength in the inte- 
rior of the Lyot plane can be found by inspection of the 
Ln,2 term in equation (|14|l . Given an aberration of the 
the form a(x 2 + y 2 ), the field strength in the interior 
of the Lyot plane is ad 2 /6, to first order. The interior 
of our coronagraph's Lyot plane fills uniformly with in- 
phase light. The effective telescope diameter describing 
the residual PSF for a circular aperture is D — 2d, so 
a fraction a 2 d 4 (l — 2d/D) 2 /36 of incident power leaks 
through the coronagraph. It is distributed in an Airy 
pattern corresponding to an aperture D — 2d in diame- 
ter, assuming an optimally undersized Lyot stop is used. 

In the absence of other aberrations, this places a 
'ghostly PSF' with a calculable intensity in the final im- 
age. The form of this PSF for a circular unobstructed 
entrance aperture is just the Airy pattern of the under- 
sized Lyot stop used in the coronagraph. This resid- 
ual PSF is not to b e confused with the spot of Arago 
l)Born fc Wold Il993l) . whose intensity varies with occul- 
ter size, depending on the relative size of an occulter in 
the pupil plane and the Fresnel zones it obscures. The 
spot of Arago's brightness fluctuates in intensity as one 
increases the pupil plane occulter size. The leaked PSF 
power we describe here does not possess undulating fluc- 
tuations, characteristic of the spot of Arago. The spot 
of Arago occurs with unaberrated wavefronts, whereas 
this leaked PSF is entirely dependent on the presence of 
aberrations in the wavefront. 



The analytical result presented here leaves no unre- 
solved questions pertaining to the numerical accuracy or 
rounding error. The approximations used in deriving this 
result become more accurate as the magnitude of the 
phase aberrations decrease, so this result is pertinent to 
coronagraphs designed for very high Strehl ratio, high 
dynamic range applications. The bandwidth of the focal 
plane mask function is the only coronagraphic parameter 
that enters into the light leak through the coronagraph 
due to inexact focus in these regimes. 

We note in passing that if there is sufficient tilt er- 
ror to merit a second order expansion of the pupil field 
strength, this second order tilt term will have a mathe- 
matically identical behavior to a small pure focus: a sim- 
ilar 'ghostly' PSF will appear in the final coronagraphic 
image. The fraction of incident power that leaks through 
the coronagraph can be calculated in a similar fashion to 
the way we estimate the power leak due to the first order 
effects of defocus. However, given the real or imaginary 
nature of the coefficients in the expansion of the elec- 
tric field in terms of the phase, the second order effect 
of tilt is in quadrature with the first order focus term 
(the former being purely imaginary, the latter real), so 
the two cannot be made to cancel each other: they will 
add power to the Lyot field (and the final coronagraphic 
PSF) in quadrature. 

If there is sufficient defocus to warrant a second order 
expansion of the pupil field, the amount of light leaking 
through the interior of the Lyot stop increases dramati- 
cally (Fig.0 second row, second column). 

6.4. Astigmatism 

Small amounts of astigmatism have no first order ef- 
fect on the band-limited coronagraph. This is easy to 
understand because of the separability of the xy polyno- 
mial, and the nature of the convolution in the definition 
of equation (|14(l . The second order astigmatic power leak 
into the Lyot plane is comparable to that of the second 
order tilt term, especially at the edges of the interior of 
the Lyot stop (Fig.0 third column). Thus, off-axis tele- 
scope designs with residual astigmatism could be used 
with Lyot stops which are undersized relative to a Lyot 
stop tailored for perfect optics. 

6.5. Coma and spherical aberration 

Leak through an unapodized Lyot coronagraph is de- 
pendent on local curvature of the wavefront, so the first 
order leak due to these terms is high, especially at the 
edges of the Lyot pupil interior. First order light leak 
into the Lyot pupil interior due to spherical aberration 
has the same form as that of the second order leak due to 
defocus (Fig. 0J, but, because of the relative quadrature 
of the first and second order field strengths, these aber- 
rations cannot be used to cancel each other by design. 

7. CONCLUSION 

Analytical studies of the mathematically simple band- 
limited coronagraph provide us with tools to help evalu- 
ate coronagraphic behavior without the need for exten- 
sive numerical investigations. The analysis also helps 
identify design strategies to reduce coronagraphic sensi- 
tivity to selected aberrations. Our results suggest that 
off-axis unobscured high dynamic range apertures can 
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be allowed more astigmatism if the Lyot pupil of a band- 
limited coronagraph is reduced in size, thereby reducing 
the demands made upon the structural stability of the 
telescope. 

The theory developed here suggests that the corona- 
graphic PSF of the nominally on-axis source can be used 
to determine best focus, by minimizing the on-axis inten- 
sity behind the focal plane mask, utilizing focus sweep 
data. If the imaging quality is sufficiently high, leaked 
power due to small tilts (i.e., small decentration of the 
target behind the focal plane mask, rather than moving 
the target out from behind the spot) can be used to de- 
termine the location of the focal plane mask. This latter 
exercise would aid astrometric calibration of the corona- 
graph. 

Further theoretical analyses of coronagraphic response 
to aberrations could assist coronagraph design and op- 
erational plans to be made, to improve the scientific 
productivity of high dynamic range space-based coron- 
agraphs dedicated to detecting and characterizing extra- 
solar planets. 
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TABLE 1 

First and second order leak vs. Strehl ratio 



99% Strehl 98% Strehl 95% Strehl 





I s * Order 


2 nd Order 


I s * Order 


2 nd Order 


1 st Order 


2 nd Order 




Power 


Power 


Power 


Power 


Power 


Power 


Tilt 





5.8 X 10" 8 





2.3 X 10" 7 





1.5 X 10" 6 


Focus 


1.8 X icr 5 


4.9 x icr 8 


3.6 X 10- 5 


2.0 X 10~ 7 


9.1 X 10~ 5 


1.3 x 10- fi 


Astig 





2.3 x icr 9 





9.4 x 10~ 9 





6.1 x 10~ 8 


Coma 


6.7 x icr 6 


5.0 x icr 9 


1.4 x 10~ 5 


2.0 x 10~ 8 


3.4 x 10~ 5 


1.3 x 10~ 7 


Sphere 


5.8 x icr 5 


3.0 x 10~ 7 


1.2 x 10~ 4 


1.2 x 10~ fi 


3.0 x 10~ 4 


7.7 x 10~ 6 



